Effect of External Noise Correlation in Optical Coherence Resonance 
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Coherence resonance occurring in semiconductor lasers with optical feedback is studied via the 
Lang-Kobayashi model with external non-white noise in the pumping current. The temporal cor- 
relation and the amplitude of the noise have a highly relevant influence in the system, leading to 
an optimal coherent response for suitable values of both the noise amplitude and correlation time. 
This phenomenon is quantitatively characterized by means of several statistical measures. 
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Despite random fluctuations usually constitute a 
source of disorder in dynamical systems, many examples 
exist in which they lead instead to an increase of order in 
the system behavior. Among these examples, stochastic 
resonance stands out for its implications in many differ- 
ent areas of science [0. In the conventional situation, 
stochastic resonance consists of an optimization due to 
noise of the response of a nonlinear system to a weak 
periodic driving [0] . But even in the absence of external 
periodic forcing, noise can be helpful in sustaining a co- 
herent oscillatory response in the system, provided the 
operation point is close to a limit cycle || or within an 
excitable regime Q. This phenomenon has been called 
coherence resonance (CR), and has been recently found 
also in bistable || and chaotic || systems. 

One of the earliest and most influential experimental 
observations of stochastic resonance was made in a laser 
system [Q. Similarly, optical systems have also provided 
in recent years clear-cut examples of excitable behavior, 
including observations in semiconductor lasers subject to 
optical feedback HH , lasers with saturable absorber jio) , 
passive nonlinear ring cavities JTl| , lasers with injected 
signal Ul, and self-pulsing lasers Following these 

studies, optical CR was predicted theoretically in the 
self-pulsing laser Jl^] and observed experimentally in the 
semiconductor laser with optical feedback M|. In this 
latter case, noise is added to the driving current of the 
laser, and gives rise to a pulsed behavior in the system, in 
the form of sudden drop-outs in the evolution of the light 
intensity. The regularity of the drop-out series initially 
increases with increasing fluctuations, and peaks for an 
optimal amount of noise. The present Communication 
is devoted to the theoretical modeling of this situation, 
making use of a rate equation system including a delay 
term, i.e. the well-know Lang-Kobayashi (LK) model 
fl6f generalized to take into account the insertion of ex- 
ternal noise into the system through the laser's pumping 
current. 

Our results show that a white-noise assumption is not 
adequate to account for the observed resonant behavior. 
In fact, such a supposition is not realistic, due to the fast 



time scales in which this system evolves (~ tens of ps), 
smaller than or on the order of the characteristic time 
scales of the fastest fluctuations that can be experimen- 
tally introduced, which are restricted by limitations of 
the electronics involved (~ GHz) . Following these con- 
siderations, we have considered a time- correlated external 
noise, and found that coherence is maximal not only for 
an optimal noise amplitude, but also for an optimal noise 
correlation time. In what follows, such a double coher- 
ence resonance is described in detail and characterized 
by suitable statistical measures. 

The LK model describes the temporal evolution of the 
slowly varying complex envelope of the electric field E(t) 
inside the laser and the excess carrier number N(t), con- 
sidering only one longitudinal mode of the solitary laser, 
and one single reflection from the external feedback mir- 
ror (i.e. multiple reflections are neglected, which is valid 
for not too large reflectivities). In dimensionless form the 
model reads pjlql: 



dE _ 1 + ia 
~dt ~ 2 



(G(E,N)-j)E(t) 



^ = 7e {cri + e(i)]iv th 



+Ke~ luJT f E(t - T f ) + y/2j3N((t) (1) 
N(t)}~G(E,N)\E(t)\ 2 , 



where 7 and j e are the inverse lifetimes of photons and 
carriers, respectively, C is the pumping rate (directly 
related to the driving current; C = 1 is the solitary- 
laser threshold), a is the linewidth enhancement factor, 
and u is the solitary-laser frequency. The last term in 
the electric-field equation represents spontaneous emis- 
sion fluctuations, with £(£) a Gaussian white noise of zero 
mean and unity intensity, and (3 measuring the internal 
noise strength. The material-gain function G(E, N) is 
given by 



G(E,N) = 



g(N(t) - N ) 
l + s\E(t)\ 2 



(2) 



where g is the differential gain coefficient and s the 
saturation coefficient. The threshold carrier number is 



1 



7V t h = j/g + No. The optical feedback is described by 
two parameters: the feedback strength k and the external 
round-trip time Tt. Finally, the external noise is repre- 
sented by the term £(t), which according to the discus- 
sion made above, is taken to be a time-correlated noise 
of the Ornstein-Uhlenbeck type, gaussianly distributed 
with zero mean and correlation 



(£(t)£(f)> = -e 



-{t-f)/r c 



(3) 



This external noise is characterized by two parameters, 
its intensity D and its correlation time r c . The variance 
of the noise is given by D/t c , and hence we will measure 
its amplitude as a = \J D /r c . 

The LK model with no external noise has been pro- 
fusely used in the past to model the dynamics of semi- 
conductor lasers subject to optical feedback. In partic- 
ular, it satisfactorily describes the appearance of drop- 
outs when the system is perturbed with electrical pulses 
over a threshold value, with the feature that the shape of 
the generated (inverted) pulses is basically independent 
of the perturbation J9|. This behavior, which is charac- 
teristic of excitable systems and appears when the laser 
is operated close to the solitary laser threshold, agrees 
with experimental observations 0. For larger pumping 
rates (but still close to the laser threshold) the inten- 
sity dropouts appear spontaneously, i.e. no external ex- 
citation is needed to produce them. Figure [l] displays 
an example of this behavior, in terms of the evolution 
of both the intensity I(t) and phase difference between 
consecutive round-trips rj(t) = <p(t) — 4>{t — r/), with 
E{t) = y/1 exp{i(f>) . The occurrence of pulses in the form 
of power drop-outs in the intensity time series can be 
clearly identified, and they are seen to correspond with 
well-defined pulses in the electric-field phase difference. 
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FIG. 1. Time evolution exhibiting drop-outs of the inten- 
sity / and the corresponding pulses in the phase difference r). 
Parameters of the LK model are: C = 1.03, 7 e = 6-l(T 4 ps-\ 
7 = 0.158 ps~\ g = 2.79 ■ 10" 9 ps" 1 , s = 3 ■ 10" 7 , a = 3.5, 
iVo = 1.51-10 8 , = 5-10" 10 ps -1 , k = 0.025 ps" 1 , Tf = 2.4 ns, 
UJTf = 2, and D — 0. 

It should be noted that the intensity time trace shown 
in Fig. l](a) has been filtered to 100 MHz, in order to 
mimic the bandwidth effect of typical photodetectors. 
For large enough filtering bandwidth (or for no filter- 
ing at all), the corresponding evolution takes the form 



of ultrashort inten sity pulses (with durations on the or- 
der of tens of ps) pM, whose envelope exhibits the low- 
frequency drop-outs shown in Fig. |l](a). 

Since the internal spontaneous-emission noise can- 
not be experimentally controlled, we now turn our at- 
tention to the effect of the external noise r/(t) upon the 
system. This effect can be understood by examining 
the mechanism behind the above-mentioned power drop- 
outs, which is well understood in the framework of the 
LK model (i.e. under the assumption of single-mode op- 
eration of the laser). This model exhibits multiple co- 
existing fixed points, which appear in pairs of solutions 
called modes and antimodes. The antimodes are saddle 
points, and most of the modes are also unstable due to a 
Hopf bifurcation |l8| . However, at least one of the modes 
(the one with maximum power) is stable. In this com- 
plex phase-space landscape, a large enough fluctuation 
may be able to take the system away from the basin of 
attraction of the stable fixed point and, upon collision 
with a neighboring antimode, produce a sudden increase 
in the phase difference [see Fig. [j](b)] which corresponds 
to a power drop-out. The corresponding escape time, 
also called activation time t a , is a random variable whose 
average decreases with the intensity of the external noise 
according to Kramers' law pj| . Following the drop-out, 
a build-up process begins in which the system undergoes 
a chaotic itinerancy around the Hopf-unstable modes, 
jumping consecutively from one to the next while being 
drifted back towards the stable maximum-gain mode [ p0[ . 
For small intensities of the external noise, the excursion 
time t e required by this process is basically independent 
of noise, and has the role of a refractory time during 
which no drop-outs can be induced. As noise intensity 
increases the escape events become more frequent, re- 
ducing the standard deviation of the interspike intervals 
accordingly. A minimum of variability occurs for an op- 
timal amount of noise when the drop-out separation is 
of the order of t e . Beyond that point, noise intensity is 
large enough to produce escapes before the build-up pro- 
cess is finished (i.e. before the stable mode is reached), 
which leads to an irregular series of pulses. This sequence 
of events is depicted in Fig. ||, which shows three time 
traces of the phase difference rj(t) for increasing ampli- 
tudes of the external noise, keeping its correlation time 
constant. In this case, the semiconductor laser is bi- 
ased at 1% above the solitary-laser threshold, a situation 
for which the system is stable in the absence of exter- 
nal noise. A small amount of noise produces infrequent 
drop-outs [Fig. ^(a)], which become more numerous and 
regular as the noise amplitude increases [Fig. @(b)]. For 
large noise strengths the pulses become increasingly ir- 
regular, both in separation and in amplitude [Fig. |2|(c)]. 
Hence, an optimal amplitude of the external noise exists 
for which the coherence of the pulsed output of the laser 
is optimal. 

In order to quantitatively characterize this effect, we 
compute the standard deviations Rg and of the nor- 
malized drop-out separation 6 = T/(T) and normalized 
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drop-out amplitude /z = A/ (A), respectively |Tq], where 
the statistical averages are performed over both time and 
different realizations of the noise. Numerically, to detect 
a drop-out we study the behavior of the phase difference 
r](t), since its variation is smoother than that of the inten- 
sity, and easier to characterize. A drop-out occurrence is 
recorded when the phase difference r\ suddenly increases 
in a fixed large amount (e.g. I2ir in our case) ||. The 
drop-out amplitude A is a measure of the total height of 
the phase pulse. 



100 



-100 
-200 
100 


-100 
-200 
100 

- 



- (a) 














\ 

1 V 'W 




!U 


V 


1 1 



03 



.(b) 

I \ (\ \ \ l\ l\ 



-200 I . i . ii . i . J i 



200 400 



600 800 

t(ns) 



1 000 1 200 1 400 



FIG. 2. Temporal behavior of the phase difference rj 
for increasing noise amplitude: (a) a = 7.36 ■ 1(T 2 , (b) 
a = 9.35 ■ 10 -2 , and (c) a = 1.60 ■ 10 _1 . In all cases we 
have kept the correlation time constant, r c = 24 ps. Other 
parameters are those of Fig. (D, except C = 1.01. 
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FIG. 3. Statistical characterization of the noise-amplitude 
coherence resonance, (a) Standard deviations Rg (full circles) 
and i? M (empty squares), and (b) joint entropy H(8,[j,) as a 
function of the external noise amplitude. The time correla- 
tion of the noise is fixed to r c = 24 ps. Other parameters are 
those of Fig. |^. 

We have computed the standard deviations Rg and 
for increasing noise amplitude, averaging up to 20,000 
drop-outs in each measure. The result is plotted in Fig. 
|§|(a), and confirms the qualitative conclusions that have 
been drawn above from Fig. ^, at least as far as the vari- 
ability of the pulse separation, Rg, is concerned. This 
quantity is a non-monotonic function of the noise am- 
plitude, being minimal for an optimal amount of noise. 
The irregularity of the drop-out amplitudes, on the other 
hand, increases monotonically with noise. This result co- 
incides with experimental observations Jl5[ , and reflects 
the fact that the frequency with which noise breaks up 
the build-up process increases steadily with the amount 
of noise added. 



In order to take into account both the drop-out sep- 
aration and amplitude simultaneously in the determi- 
nation of the signal's regularity, it is useful to define 
a joint entropy H(8,fi) of the two quantities, where 
H = — ^PlogP, with P the joint probability density 
of the two random variables pi| . For the case of two 
Gaussian independent random variables the following re- 
lation holds [pf: exp(H(9,fi)) = 2ireRgR^. We will as- 
sume that this result is approximately valid in our case, 
and compute the joint entropy accordingly. The result 
is given in Fig. ^|(b), which shows again a maximum 
regularity of the drop-out series for an optimal noise am- 
plitude, this time taking into account both the pulse sep- 
aration and amplitude. In fact, the minimum is in this 
case more clearly defined. 

We note that all results presented so far have been com- 
puted for a fixed correlation time of the noise r c = 24 ps, 
on the order of the fast time scale of the deterministic 
dynamics. In fact, as the white-noise limit is approached 
the amount of noise necessary to obtain similar effects 
climbs up to unreasonably high values. The reason is 
that the carrier dynamics acts as a frequency filter for the 
external noise [see equation for N(t) in (Q)], which also 
prevents the system from responding to high-frequency 
modulations of the pump current. Therefore, most of 
the power of a white noise has no effect upon the sys- 
tem dynamics, and the noise intensity needs to be very 
large in order to have a noticeable influence (a similar 
effect has been observed in periodic-modulation studies 
p2| ). In the opposite frequency limit a similar situation 
occurs: for low-frequency forcing the carrier dynamics 
has enough time to follow the modulation, and the sys- 
tem responds simply with a modulated output. Only for 
intermediate frequencies will the external forcing be able 
to influence the drop-out statistics and enhance the co- 
herent response of the system. In order to verify this 
conjecture, we now fix the amplitude a = \/D/t c of the 
external noise and analyze the behavior of the system for 
an increasing correlation time of the Ornstein-Uhlenbeck 
noise defined by Eq. (||). The result is shown in Fig. 
H for three different values of r c . It can be seen that 
the regularity of the pulsed time series is maximal for 
intermediate values of the noise correlation time. 



We quantify again the qualitative observation made in 
the previous paragraph by computing the standard devi- 
ations Rg and R^, and the joint entropy H(6,fi) of the 
normalized drop-out separation and amplitude. The re- 
sults, shown in Fig. (|^), exhibit the same behavior as in 
the case of an increasing noise amplitude. Coherence of 
the pulsed behavior is in this case maximal for a corre- 
lation time of the noise r c ~ 30 ps. This behavior can 
be interpreted as a resonance with the fast deterministic 
dynamics of the system. A similar resonance has been 
recently observed experimentally in a chemical excitable 
medium |G3|. 
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FIG. 4. Temporal behavior of the phase difference r\ 
for increasing noise correlation time: (a) r c = 15.8 ps, (b) 
t c = 57.6 ps, and (c) r c = 153.2 ps. In all cases we have kept 
the noise amplitude constant, a — 0.079. Other parameters 
are those of Fig. B. 
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Statistical characterization of the 
noise-correlation coherence resonance, (a) Standard devia- 
tions Re (full circles) and (empty squares), and (b) joint 
entropy H(6, fx) as a function of the noise correlation time. 
The noise amplitude is fixed to a = 0.079. Other parameters 
are those of Fig. ^. 

In conclusion, we have shown the existence of dou- 
ble coherence resonance in an excitable optical system 
(a semiconductor laser with optical feedback) driven by 
an external correlated noise. The assumption of a non- 
delta temporal correlation of the noise is fully meaning- 
ful, given the fast characteristic time scale of the deter- 
ministic dynamics of the system. The pulsed response 
of the laser, which takes the form of intensity drop-outs, 
exhibits a maximal coherence for optimal values of both 
the amplitude and correlation time of the external noise. 
These results agree satisfactorily with previous experi- 
mental observations. Other recent investigations have 
analyzed the influence of correlated noise in model sys- 
tems exhibiting coherence resonance 24 , but in those 



cases the corresponding deterministic system lacked a 
natural second time scale (different from the optimal 
pulse separation) which could account for the resonance 
with respect to the noise correlation time. In the present 
case, however, such a time scale does exist naturally in 
the system. 
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